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We investigate the evolution of non-adiabatic collapse of a shear-free spherically symmetric stellar 
configuration with anisotropic stresses accompanied with radial heat flux. The collapse begins from 
a curvature singularity with infinite mass and size on an inhomogeneous space-time background. The 
collapse is found to proceed without formation of an even horizon to singularity when the collapsing 
configuration radiates all its mass energy. The impact of inhomogeneity on various parameters of 
the collapsing stellar configuration is examined in some specific space-time backgrounds. 
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I. INTRODUCTION 



One of the most fundamental problems in general relativity is the construction of realistic models describing various 
evolutionary stages of a star collapsing under its own gravity. When a star exhausts all its thermonuclear fuel, it 
, cannot withstand the gravitational pull and starts collapsing. As it contracts appreciably, the density increases and 
Q H ' at a sufficiently high density it produces non-thermal pressure via degenerate fermions and particle interactions to 
support it against further collapse and becomes a 'compact star'. Stable stellar configurations like neutron stars are 
(^JT) 1 the end products of such radiative collapse processes. A massive star, however, can not come to a stable stage by 
l— ~~ l \ such processes. In the absence of any mechanism that can withstand the gravitational pull, the general relativistic 
prediction is that such a collapse must terminate into a space-time singularity. In view of the Cosmic Censorship 
Conjecture(CCC), such a singularity must be covered within its event horizon of gravity [lj. However, there are several 
■ counter examples where a naked singularity is more likely to be formed. In fact, there is yet no established theory 
t-H ' available governing the formation of either a black hole or a naked singularity (for a recent review see |2| and references 
therein). The physics of evolving dynamical systems, therefore, continues to generate a great deal of interest in various 
fields of astrophysics and cosmology even today. 
\Q \ To understand the nature of collapse of a self gravitating body, one needs to provide an accurate description of the 
exterior and interior space-times of the collapsing body. It is also necessary to find out the appropriate boundary 
conditions joining the two regions. The theoretical understanding of gravitational collapse was first initiated by 
Oppenheimer and SnyderQ who considered the contraction of a highly idealized spherically symmetric dust cloud. 
The exterior space-time of the collapsing dust cloud was described by the Schwarzschild metric and the interior 
space-time was represented by a Friedman-like solution. Later on, the VaidyaQ metric corresponding to the exterior 
gravitational field of a stellar body describing an outgoing null fluid gave a tremendous impetus in this direction. 
Making use of the Vaidya[H solution, SantosQ proposed a procedure to obtain the description of the interior space- 
time of a spherically symmetric radially shrinking distribution of non-adiabatic fluid. Several stellar models (see 
for example, [6 34] and references therein) have been obtained and examined critically using Santos's approach. 
The procedure has been found to be useful to examine the impact of various factors such as shear, inhomogeneity, 
anisotropy, electromagnetic field and various dissipative processes on the evolution. In the absence of any established 
theory governing the nature of collapse, such investigations have been found to be very useful to understand the 
dynamic behaviour of a gravitationally collapsing system. 

The objective of the present work is to formulate a framework to study the nature of collapse of an inhomogeneously 
distributed anisotropic source. The dynamical behaviour of a gravitating system is expected to be influenced by density 
inhomogeneity and anisotropy. Amongst many other issues, these factors have been found to play crucial role in our 
theoretical understanding of the nature of singularity and its genericity (see [l], d, [3oT - l39j and references therein). 
Eardley and Smarr 40] have shown that inhomogeneous distributions may lead to a naked singularity in contrast 
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to a homogeneous distribution where a black hole is more likely to be formed. For a spherically symmetric dust 
cloud, Mena et aZl4ll have examined the role of inhomogeneity and anisotropy in Lamaitre-Tolman-Bondi collapse. 
Hererra et a/[l6|, EH have examined the role of various factors contributing to inhomogeneity of matter distribution, 
its evolution and development of anisotropy. It has been observed that tidal forces tend to make a gravitating system 
more inhomogeneous which may lead to the formation of a naked singularity while for a homogeneously distributed 
body it is more likely to form a black hole. Similar observations may be found in [2lj ] where it has been shown that 
if an initially static star undergoing collapse has a homogeneous distribution it ends up with the formation of a black 
hole. However, if an anisotropic star contains an inhomogeneous distribution of matter before the collapse sets in, 
then for a shearing model of collapse, the black hole is never formed as the condition for the formation of apparent 
horizon is never satisfied under such a condition (2l|. 

Pressure anisotropy also plays a significant role in the construction of a dynamically evolving system. Elaborate 



discussions on the microscopic origin of anisotropic stresses in stellar bodies may be found in 14(,[43(. Amongst 
many other factors, electromagnetic field and shear have been found to contribute significantly to the generation of 
anisotropic stresses. The effect of charge on non-static models of gravitationally bound objects have been investigated 
by many authors @, Gil El EI EMI SS SI IH IH 0, Hi] and references therein) . Shear is also a generator 
of anisotropic stresses [28|. Effects of shear and anisotropy in the dynamical behaviour of a self gravitating system 
have been analyzed by Govinder et al[46] . It has been observed that if a star has isotropic pressure before the collapse 
sets in, anisotropy may develop at a later stage due to the presence of shear (20l [22l|. 

In our work, we have developed a model describing a spherically symmetric inhomogeneous anisotropic fluid config- 
uration radiating away its energy in the form of radial heat flux and shrinking in size as the collapse proceeds. In our 
construction, we have assumed that the collapsing star is charge neutral and the back ground space-time is shear-free. 
However, the energy-momentum tensor corresponding to the fluid distribution filling the interior of the collapsing 
star has been assumed to be anisotropic, in general. The model developed here is a generalization of an earlier model 
presented by Banerjee et a/[48| describing the collapse of a homogeneously distributed isotropic fluid. Schafer and 
Goenner[4!| used the Banerjee et a/[48[ model to investigate the collapse of a homogeneously distributed isotropic 
fluid configuration and also critically a naly zed various features of the model parameters involving the background 
space-time. In the Schafer and Goenner[49| model, the collapse begins at time t = — oo with both infinite mass and 
radius and contracts to a point at time t = without forming an event horizon. The horizon is never formed in 
this set up because the rate of collapse is counter balanced by the rate at which energy is dissipated to the exterior 
space-time of the star. 

In our construction, we have traced the evolution of the collapse on the background of space-time obtained by 
introducing an inhomogeneous perturbation in the Robertson- Walker space-time. Inhomogeneity in the geometry of 
the background space-time may be interpreted in the following way. We have considered the 3-sub-space of the 4-D 
manifold as having a geometry of a 3-spheroid rather than a 3-sphere as in 48]. We have examined the impacts of 
inhomogeneous nature of background space-time and anisotropic stresses on the collapse by comparing the behaviour of 
physical parameters in our set up to the behaviour of corresponding parameters in Banerjee et a/[48] model, admissible 
as a special class in our model. We have also explored the thermal behaviour at the interior of the collapsing star 
by considering two different space-time backgrounds. We have not incorporated bulk viscosity and electromagnetic 
field in our model as these factors are likely to be absorbed into the anisotropic stresses. We have also considered 
a shear-free model so as to keep the governing field equations reasonably simple and tractable. Moreover, since the 
Schafer and Goenner[4^| model was assumed to be shear- free, this simplifying assumption helps us to examine the 
impact of inhomogeneity directly by comparing our results to its homogeneous counterpart. 

The paper has been organized as follows: In Section 2, we have presented an inhomogeneous generalization of the 
Banerjee et a^48| model by introducing a perturbation in the background space-time for a spherically symmetric 
anisotropic fluid undergoing non-adiabatic radiative collapse. Stipulating the boundary conditions across the surface 
separating the stellar configuration from the Vaidya[4j space-time describing its exterior filled with outgoing radiation, 
we have solved the surface equation which governs the evolution of the collapse. In Section 3, we have examined 
the implications of inhomogeneity on the collapse by comparing the evolutions of the physical quantities to the 
homogeneous model discussed in [49( . Bounds on the model paramters based on physical requirements and evolution 
of temperature have been analyzed in this Section. Finally, some concluding remarks have been made in Section 4. 
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II. INTERIOR SPACE-TIME 



The conformally flat space-time metric formulated by Maiti[47 
and rotation-free fluid with heat flux as source, has the form 



representing a spherically symmetric shear-free 



ds 2 
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The line element (HJ reduces to Robertson- Walker metric for a(t) — k — 0. Making use of the metric (JTJ, Banerjee et 
q/[48| presented a simple model for a collapsing body with outgoing radiation by assuming a(=^ 0) and k as remaining 
constants during collapse. The collapsing matter in this set up was assumed to be a homogeneously distributed 
isotropic fluid. Various aspects of this model have been extensively examined by Schafer and Goenner[49( putting 
constraints on model parameters complying with various physical plausibility requirements. Following Schafer and 
Goenner[49|. we express the space-time metric of Banerjee et a/[48] model in standard coordinates as 



dsi 



-(C- ^l-kr 2 ) 2 dt 2 +R 2 (t) 



dr 2 



kr 2 



(2) 



where, C and k(^ 0) are constants. 

To generalize the Banerjee et aZ[48] model, we assume that the collapsing configuration is comprised of an inhomo- 
geneous distribution of anisotropic fluid with its background space-time having the form 



ds 2 



-A 2 {r)dt 2 + R 2 {t) 
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(3) 



Here, the metric potential Ao(r) and the scale factor R(t) are undetermined metric functions and A is a parameter 
measuring departure from homogeneous geometry. The t = constant hyper surface of the space-time (|3]) has a 
geometry of a ■''> spheroid ">() "> 1 representing a perturbation from that of a 3-sphere. The energy-momentum tensor 
of the fluid with anisotropy in pressure filling the interior of the collapsing body is written explicitly in the form 



T af 3 = {p+Pt)u a Up +Pt9a/3 + (Pr -Pt)XaXp + QaUfj + qpU a . 



(4) 



Here p represents the energy density, p r is the radial pressure, pt is the tangential pressure, \ a is a un it space like 
four vector along the radial direction, u a is the 4-velocity of the fluid and q a = qS" is the heat flux vector which 
is orthogonal to the velocity vector so that q a u a = 0. Einstein's field equations in view of © and (HJ) lead to the 
following system of four independent equations 
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where, 
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We have used in above the system of units rendering G = c = 1. Combining Eqs. ©-(O, a time-independent 
differential equation of the form 
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may be obtained if it is assumed that the anisotropy evolves as 



8Tr{pt-p r ) = A(r,t) 



R 2 (ty 



By making a transformation x 2 = 1 — kr 2 and using the expression for Bo(r), Eq. (|10|) may be rewritten as 
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In terms of a new dependent variable 



*(a?) = A (x)[l + Afc - Afca: 2 ] -1 / 4 , 



Eq. I|12p assumes the form 
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^/ can be determined from Eq. (|14[) if the nature of anisotropic parameter 5(r) is known. We shall examine the effect 
of inhomogeneity on the evolution of the collapse using the simple tractable solution 



* = C + Dx, 



of Eq. ([T4"|) , obtained by stipulating 
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Note that the anisotropy parameter 5(r) is regular at all interior points of the configuration. In Eq. (|T5j) . C and D 
are arbitrary constants of integration. Combining Eqs. (|13[) and (|15p . the metric function A (r) is obtained as 



A (r) = (1 + Afc 2 r 2 ) 1/4 (C + L> ^1 - kr 2 ). (17) 
Consequently the space-time of the collapsing shear-free stellar body with anisotropic stresses is described by metric 

~ 1 + Afcr 2 



-(1 + Afc 2 r 2 ) 1 / 2 (C + D^Jl - kr 2 fdt 2 + R 2 (t) 
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(18) 



The Banerjee et aZ[4a| model is a sub-class of (fT8|) which follows on setting A = and D = — lj49|. We shall examine 
the impact of the inhomogeneous nature of the background geometry due to presence of A without any loss of generality 
by stipulating D = — 1 . 



A. Determination of R(t) 

The evolution of the collapse in the stellar body with interior space-time metric (fT8)) is governed by the function 
R(t) which has role of a scale factor in the process. It is determined by the boundary conditions across the boundary 
surface S of the configuration as it shrinks in size. The space-time outside the collapsing body will be appropriately 
described by the Vaidya[iJ metric 

ds\ = - (l - dv i _ 2dvdf + r 2 dn 2 . (19) 

Here v denotes retarded time and m(v) represents the total mass of the collapsing star which is expected to be a 
function of the retarded time v. Following the method presented by Santos[5j, we write the matching conditions 
linking smoothly the interior and the exterior space-times across the boundary 3-space-time of the evolving system as 
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From (|20p it follows that the mass of matter enclosed within the boundary surface r < rs will be 
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Eqs. ((BJ, (jH) and the boundary condition (|2T|) at r = rs lead to 
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A simple solution of Eq. ([23]) will then be [4' 



where, 



n = a ± Va 2 - 2/3. 

For collapsing configurations we must have R < 0. Hence n should be negative implying 

n = a — \/ a 2 — 2/3. 
The constants C, fc and A should be chosen to ensure a 2 > 2(3 and n < 0. 
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III. PHYSICAL ANALYSIS 



Einstein's field equations lead to the following explicit expression for the dynamical variables of matter density, the 
radial and transverse pressures and the heat flux parameter associated with the collapsing stellar structure 
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These physical parameters have zero values at t — > — oo; the density and the two pressures evolve with time as ~ 1/t 2 
while the heat flux evolves as ~ 1/t 3 . The proper radius r p = [R(t)r]z = ntr^, is infinite when collapse begins, positive 
at any later instant t and shrinks to zero at t = 0. The mass at any instant t within the boundary radius has 
expression 



m(v) = m{rz,t) = — — 



(1 + A)fcr| n 2 r| 

(1 + Xkr 2 ,) + (C - V /T _ fc4) 2 (l + Afc 2 r|jV 2 



(33) 



which shows that the collapse begins with an infinite mass and size of the configuration at t — > — oo and it evaporates 
completely as the epoch t = approaches. As in [56J, the collapsing matter radiates all the mass energy when the 
singularity is reached as proper boundary radius of the configuration r p — [rR(t)]^ = at t = 0. 
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A. Nature of singularity 



To understand the nature of singularity, in general, it is necessary to verify the existence of radial null geodesies 
emerging from the singularity. The existence of a naked singularity is characterized by the absence of apparent horizon 
as the collapse approaches the central singularity r = 0. However, such a test is not required in our model since the 
ratio of (2m(r,t)/rR)s is independent of time. One can always suitably choose the model parameters so that the 
ratio (2m(r, t)/rR)s < 1. The function 1 — [2m(r, t)/(rR)]z will then remain positive throughout the collapse process 
which indicates that the collapse would terminate into a naked singularity. The Ricci curvature, in our model, turns 
out to be 

_ 6r 2 n 2 B 2 - 2A 2 (B - 2rB' ) + 2rA [A' B' - B (2A> + rA'>)} 

r 2 A 2 B 2 n 2 t 2 ' { ' 

which show that the curvature diverges as 1/t 2 . This suggests that it is a weak curvature singularity. 



B. Energy conditions 

A physically reasonable solution should satisfy certain energy conditions, namely, (1) the null energy condition 
(p > 0); (2) the weak energy condition (p—p r > 0, p— pt > 0); and (3) the dominant energy condition (p—p r — 2p t > 0). 

From Eq. (|29|) - (f5T]) . it is easy to show that p, p r , p t > and p' , p' r , p' t < 0, if the conditions A > and k > 
are satisfied simultaneously. Eq. (|30[) shows that p > p r . Note that at r — 0, p r — pt and both p r and pt decrease 
radially outward. We, therefore, may conclude that the null and weak energy conditions are satisfied in our model. 
It is difficult to verify the dominant energy condition unless numerical techniques are invoked. However, we note that 
by setting R(t) = 1 in our model, it is possible obtain an initial static configuration. The t — constant hypersurface 
of the initial static configuration describes a space-time geometry which is spheroidal in nature. The corresponding 
static configurations have been found to comply with the energy conditions mentioned above (see [55| and references 
therein). Since p, p r and pt evolve as 1/t , it is expected that the fulfillment of the energy conditions will be sustained 
throughout the collapse. 



C. Bounds on model parameters 



The rate of expansion, in our model, has the explicit expression 

Q = u a a = = (35) 

1 A R t(l + Afc 2 r 2 ) 1 /4(C- Vl-fcr 2 ) V ' 

and should be negative since our model describes a contracting body. Therefore, during the whole collapse process 
from time t = — oo to t = 0, we must have 



(1 + AfcV) 1/4 (C - Vl - kr 2 ) > 0. (36) 

Eq. (|3"6")l implies a lower bound on C such that C > 1. Eq. (fBl)]) also follows from the requirement Ao(r) > 0. The 
constraints on the model parameters obtained in view of physical viability are specified below: 

(a) From Eq. ([29]) . the central density is obtained as 

3fc(A + l) 3 



n 2 i 2 i 2 (C-l) 2 ' 

Therefore, for a positive energy density, we must have k > 0. Moreover, the requirements p, p r , pt > and 
p', p' rl p' t < 0, demand that A > 0. 



(b) For < r < rs, the factor vl — kr 2 appearing in the metric function Aq(t) will be real if < k < 1/rf;. 

(c) As discussed earlier, the lower bound on C is given by C > 1. The upper bound on C may be obtained by noting 

that a 2 > 2(3 which ensures that n in Eq. (|28|l remains real. Thus, bounds on C may be expressed as 



1< C < 



8y/l - fer| + 4i/l + Afcr 2 (l + Afc 2 r|) 2 + - kr- 



4 + 6A 
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where. 

a = 4-4/c + 22/s 2 r|+8Afc 2 r|-5Afc 3 r| + 15Afc 4 r|+4A 2 fc 4 r|, 
b = 4-4fc+12A: 2 r|+8Afc 2 r|-5Afc 3 r| + 9AA: 4 4+4A 2 fc 4 r 4 ] . 

(d) For collapsing configurations we must have R < which implies that n < 0. 

Though, more stringent constraints on the model parameters can not be ruled out, the constraint equations obtained 
previously by Schafer and Goenner|49| may be regained by stipulating A = in this model. The parameter A will 
take on value such that all the above mentioned conditions are satisfied simultaneously. In Fig. ([T]), we have shown 
the range of values for k and C consistent with the constraints which can ensure regular behaviour of the physical 
parameters. Note that the condition (d) yields a more stringent bound on the parameters as compared to the condition 

We have checked the behaviour of the physical parameters like energy density, radial and tangential pressures by 
assuming a set of values of A, k and C consistent with the above constraints. Fig. ©-([T]) describe the behaviour of 
various physical parameters and indicate that they all comply with the requirement of the regularity. In the figures, 
we have assumed that the collapsing object has an initial radius rs(t — > — oo) = r s = 1 as in [4!|. The proper 
radius r p = R(t)rs = ntrs, shrinks to zero at t — as R(t) = at t = 0. Therefore, a choice of r s is possible for 
convenience and will not have any impact on the nature of collapse. The choice r s = 1 has the effect of changing 
the length measuring scale only. The bounds on the model parameters are valid for all < rs < oo including 
rs{t -> -oo) = r s = 1. 

In Fig. ©-(H]), we have shown the evolution of the central values of the energy density, radial and tangential 
pressures, respectively, by assuming a set of values consistent with the constraints (we have assumed A = 5, k — 0.2, 
C = 1.2) for an inhomogeneous anisotropic distribution. Setting A = 0, k = 0.2, C = 2.5, behaviour of the 
corresponding parameters for a homogeneous isotropic distribution have been shown in Fig. l[5 ]) -([7 ]) . The plots clearly 
indicate that the physical parameters start with zero value at t — > — oo and increase till the singularity is reached 
at t = and = 0. Since anisotropy vanishes at the centre, evolution of the radial and tangential pressures show 
identical behaviour though explicit expressions for p r and p t in Eqs. (f50"|) - (pTlT) clearly indicate their inequality for 
r =/= 0. If A = 0, p r = pt at all interior points of the star which implies that pressure isotropy condition can be regained 
by setting A = 0. In Fig. |(5J) -([9]), we have shown the variations of the total mass m(r^,t) for A = and A = 5, 
respectively. Evolution of the scale factor R(t) has been shown in Fig. (fTU|) for two different values of A(= 0, 5). The 
A = case corresponds to the Schafer and Goenner[4^| model describing the collapse of a homogeneous fluid model 
with isotropic fluid pressure. 



D. Thermal behaviour 



Let us now investigate the evolution of temperature of the collapsing star. In view of extended irreversible ther- 
modynamics, the relativistic Maxwell-Cattaneo relation for temperature governing the heat transport within the 
collapsing matter in the truncated Israel- Stewart theory [57rf6Cl | has the form 

r(g a P + u a uP)u s q^ s + q a = -k( 9 ^ + u a u^)[T p + Tup], (37) 

where k(> 0) is the thermal conductivity and r(> 0) is the relaxation time. In view of the line element (ITBl . Eq. (|3"7} 
reduces to 

T^-(qRB ) + qRA Q Bo = -k-^--^-(A T) (38) 
at Rr>o dr 



The relativistic Fourier heat transport equation may be obtained by setting r = in ([38]). To get a simple estimate 
of the temperature evolution, we set r = in our calculations. For r = 0, combining Eqs. © and (13"8"1) . we obtain 

8 ™(^° T )' = ^f> ( 39 ) 

Following an earlier treatment [60j. we choose the thermal conductivity parameter in the form k = jT u , where 7 and 
uj are constants. Eq. (l39l) . then, yields 

OA' rp-U 

SnjiAoTY^-V—, (40) 
An t 
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where we have used Eq. (|26|) . Integrating the above equation, we get 

8^A T^ = (41) 
The surface temperature at any instant for the special case (u> = 0) in our model thus turns out to be 



4^(1 + Afc 2 r|) 1 /4(C - Dy/1 - fcr 2 ) ' 

Obviously, the surface temperature has zero value at t — > — oo and it evolves as 1/t. Time evolution of the surface 
temperature for inhomogeneous (A = 5) and homogeneous (A = 0) distributions has been shown in Fig. (|11[) . An 
interesting feature of our model is that the surface temperature for the inhomogeneous, anisotropic distribution 
appears to be less than that of its homogeneous, isotropic counterpart. However, evolution of the temperature shows 
identical behaviour in both the cases. 



IV. DISCUSSION 



We have generalized the Banerjee et aZ[48l] model and examined the impacts of inhomogeneous nature of background 
space-time and anisotropic stresses on the collapse by comparing the behaviour of physical parameters in our set up to 
the behaviour of corresponding parameters of a collapsing homogeneous, isotropic stellar configuration of Schafer and 
Goenner[4{| model, admissible as a particular class. We have considered two cases setting the parameter A = 0, 5, 
for studying the evolution of the collapse using numerical procedures. The parameter A is a geometrical parameter 
which is a measure of inhomogeneity. We have derived new bounds on the model parameters and the bounds on the 
model parameters found in [49( may be regained simply by setting A = 0. 

Though, in general, the physical parameters evolve differently for two different back ground space-times, an in- 
homogeneous perturbation of the background space-time of the collapsing stellar configuration seems to have very 
little impact on the gross features of the evolving system. Apparently departure from homogeneity does not have any 
drastic impact on the evolution of collapse except for a change in the collapse rate by a factor (1 + Afc 2 ^ 2 ) 1 / 4 . The 
possibilities of an in-depth analysis in different set up giving conflicting results cannot be ruled out. 
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FIG. 1: Range for the parameters k and C for A = (blue shaded region which includes the green shaded region) and A = 5 
(green shaded region) obtained by imposing the condition n < 0. For A = and 5, the blue and green curves respectively give 
the bounds on the model parameters obtained from the condition (c). 
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FIG. 2: Evolution of the energy density (87rp) at the centre (r = 0) for an inhomogeneous and anisotropic distribution (A = 5, 
k = 0.2, C — 1.2). 




FIG. 4: Evolution of the transverse pressure (8npt) at the centre (r = 0) for an inhomogeneous and anisotropic distribution 
(A = 5, k = 0.2, C = 1.2). 
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FIG. 5: Evolution of the energy density (87rp) at the centre (r = 0) for a homogeneous and isotropic distribution (A = 0, 
fc = 0.2, (7 = 2.5). 
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FIG. 6: Evolution of the radial pressure (87rp r ) at the centre (r = 0) for a homogeneous and isotropic distribution (A = 0, 
fc = 0.2, (7 = 2.5). 
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FIG. 7: Evolution of the transverse pressure (8npt) at the centre (r = 0) for a homogeneous and isotropic distribution (A = 0, 
fe = 0.2, (7 = 2.5). 




FIG. 9: Evolution of the mass function m(rs,f) for an inhomogeneous and anisotropic distribution (A = 5, k — 0.2, C = 
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FIG. 10: Evolution of Scale factor R(t) for homogeneous (A = 0) and inhomogeneous (A 7^ 0) distributions. 




